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0^\ Abstract 

^Q \ The EPDifF equation (or dispersionlcss Camassa-Holm equation in ID) is a well known ex- 

ample of geodesic motion on the Diff group of smooth invcrtible maps (diffeomorphisms) . Its 

pg \ recent two-component extension governs geodesic motion on the semidirect product Diff (S) ^, 

where T denotes the space of scalar functions. This paper generalizes the second construction 
to consider geodesic motion on Diff (S)g, where g denotes the space of scalar functions that 

jy^ • take values on a certain Lie algebra (for example, g = J^ g) so(3)). Measure-valued delta- 

pH I like solutions are shown to be momentum maps possessing a dual pair structure, thereby 

extending previous results for the EPDiff equation. The collective Hamiltonians are shown 
to fit into the Kaluza-Klein theory of particles in a Yang-Mills field and these formulations 
are shown to apply also at the continuum PDE level. In the continuum description, the 

-^ , Kaluza-Klein approach produces the Kelvin circulation theorem. 
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1 Introduction 

1.1 Singular solutions in continuum mechanics 

Singular measure-valued solutions arise in the study of many continuum systems. A famous 
example of singular solutions in ideal fluids is the point vortex solution for the Euler vorticity 
equation on the plane. Point vortices are delta-like solutions that follow a multi-particle dy- 
namics. In three dimensions one extends this concept to vortex filaments or vortex sheets, for 
which the vorticity is supported on a lower dimensional submanifold (ID or 2D respectively) 
of the Euclidean space M'^. These solutions form an invariant manifold. However, they are not 
expected to be created by fluid motion from smooth initial conditions. Singular solutions also 
exist in plasma physics as magnetic vortex lines and in kinetic theory as single-particle solutions, 
see, e.g. [H] for discussion and references. 

Another well-known fluid model admitting singular solutions is the Camassa-Holm (CH) 
equation, which is an integrable Hamiltonian fluid equation describing shallow-water solitons 
moving in one dimension ^. In the dispersionless limit, the CH solitons are singular solutions 
whose velocity profile has a sharp peak at which the derivative is discontinuous. These are the 
peakon solutions of the CH equation. In higher dimensions the dispersionless CH equation is often 
called EPDiff , which is short for Euler- Poincare equation on the diffeomorphisms |181 [T9] . EPDiff 
also has applications in other areas such as turbulence [8] and imaging [21^ [22] . EPDiff and the 
Euler fluid equation share the important property of being geodesic motions on diffeomorphism 
groups. In particular, EPDiff involves the entire Diff group, while its restriction to the volume- 
preserving transformations Diffyoi yields the Euler fluid equation. Moreover, the EPDiff equation 
has the additional interesting feature of showing spontaneous emergence of singular solutions 
from any confined initial velocity configuration. In ID this result follows from the steepening 
lemma [1]. The dynamical variable is the fluid velocity (or momentum) and the peaks of the 
singular solutions follow the Lagrangian trajectories of fluid particles. Recent work has proven 
that the singular solutions of EPDiff represent a certain kind of momentum map |18j . This 
result answers the fundamental question concerning the geometric nature of these solutions. In 
particular, they are equivariant momentum maps admitting a dual pair representation similar 
to that found for vortex filaments in |29j . 

The Camassa-Holm equation in ID also possesses a two-component integrable extension 
(CH2) [6l[71[23], which represents geodesic motion on the semidirect product Diff(s).F, where 
J^ denotes the space of real scalar functions. This system of equations involves both fluid 
density and momentum, and it also possesses singular solutions in the latter variable. CH 
peakon solutions form an invariant subspace of CH2 solutions for the case that the density 
vanishes identically. A change in the metric allows delta-like singularities in both variables, not 
just the fluid momentum, although such a change may also destroy integrability [22]. This 
paper investigates the dynamics and the geometric nature of singular solutions of the geodesic 
equations on semidirect product Lie groups. In particular, it shows that these solutions are still 
momentum maps; so that all the results for momentum maps in the case of EPDiff also hold for 
this system. 

1.2 The Camassa-Holm and EPDiff equations 

The dispersionless Camassa-Holm equation is a geodesic flow on the infinite-dimensional Lie 
group of smooth invertible maps (diffeomorphisms) of either the real line or a periodic interval. 
This geodesic flow exhibits spontaneous emergence of singularities from any confined smooth 
initial velocity profile. The CH equation is a 1-1-1 partial differential equation (PDE) for the ID 



fluid velocity vector u as a function of position x G M and time t G M and is written as [1] 

Ut + 2kUx - Uxxt + 3uUx = 2UxUxx + UUxxx ■ (1-1) 

The present work focuses on the CH equation (jl.ip in the dispersionless case for which k = and 
considers periodic boundary conditions or sufficiently rapid decay at infinity, so that boundary 
terms do not contribute in integrations by parts. Being geodesic, the CH equation can be 
recovered from Hamilton's principle with a quadratic Lagrangian, 

5 L{u)dt = with L{u) = - I u{x) {I - dl) u{x) dx . 

In particular, it follows from the Euler-Poincare variational principle defined on X(M) = TgDiff (R), 
the Lie algebra of the Diff group, consisting of vector fields on the line [19]. The corresponding 
Lie-Poisson Hamiltonian formulation follows from the Legendre transform 

r r 

m = -^ = u — Uxx ^ u = (1 — d^)~ m , 
5u 

in which m € X*(M). Here X*(M) is the space of one- form densities, which are dual to the vector 
fields on the line under the L^ pairing. After the Legendre transformation, the Hamiltonian 
becomes 

Him) = - I mil — drr)" mdx =^ u = - — 
2 7 5m 

and the corresponding Lie-Poisson form of equation (jl.ip with k = is ^ 

mt + unix + 2uxm, = . 

The main result for this equation is its complete integrability, which is guaranteed by its bi- 
Hamiltonian structure [1]. Another important feature, called the steepening lemma [3], provides 
the mechanism for the spontaneous emergence of the singular solutions (the peaked solitons, or 
peakons mentioned earlier) from any confined initial velocity distribution. 

EPDiff. Except for integrability, these one-dimensional results can be generalized to two or 
three dimensions, in which the equation becomes EPDiff, namely, 

dtm — u X curlm + V(u • m) -|- m(divu) = . (1-2) 

The EPDiff Hamiltonian on X*(M^) arising from this generalization is given by 

1 /" 1 

H{m) = - / m • (1 - a^A)"^mdx =: -||mf , (1.3) 

in which a is the length-scale over which the velocity is smoothed relative to the momentum via 
the relation 

u = (l-a^A)"^!!!. (1.4) 

The singular solutions in higher dimensions are written in the momentum representation as 

N 

m(x, t) = J2 Pi(^. *) '5(x - Qi(^' 0) ds , (1.5) 

1=1 -^ 

where s is a variable of dimension k < 3. These solutions represent moving filaments or sheets, 
when s has dimension 1 or 2, respectively. 



A further generalization replaces the kernel that defines the norm of ||m|| in p.3|) via the 
relation (|1.4p with the general convolution 

G * m= f G{x-x) m(x') dx , (1.6) 

involving an arbitrary Green's function, or kernel, G. This generalization produces the metric, 

i?(m) = iym-(G*m)dx=:i||m||2,, (1.7) 



which is the Lie-Poisson Hamiltonian for EPDiff (jl.2p and is, thus, invariant under its evolution. 
The dynamics of (Qj, Pj) with i = 1, . . . ,N for the singular pulson solutions in (II. Sp is given by 
canonical Hamiltonian dynamics with the Hamiltonian 

obtained by evaluating the Lie-Poisson Hamiltonian for EPDiff on its singular solution set. 

Definition 1.1 (Momentum map [26]) Given a Poisson manifold (i.e. a manifold P with a 
Poisson bracket {■,■} defined on the functions J-{P)) and a Lie group G acting on it by Poisson 
maps (that is, via the Poisson structure), a momentum map J : P — > g* «s defined by the 
following formula, 

{F{p),{J{p),0}=^p[F{p)] VFE^(P), V^Gfl. 

Here J^{P) denotes the functions on P, g is the Lie algebra of G and ^p is the vector field given 
by the infinitesimal generator of the action of G on the manifold P by Poisson maps. 

Theorem 1.2 The singular solution \1.5(j is a momentum map IIS^ . 

To illustrate this theorem, fix a /c-dimensional manifold S immersed in M" and consider 
the embedding Qj : 5 — > M". Such embeddings form a smooth manifold Emb(S', M") and thus 
one can consider its cotangent bundle (Qj,Pj) G T*Emb(5', M"). Consider Diff(M") acting 
on Emb(5, M") on the left by composition of functions (sQ = 5 o Q) and lift this action to 
r*Emb(S', M"). This procedure constructs the singular solution momentum map for EPDiff, 

J :r*Emb(5,M")^X*(M") with J{Q,P) = f F{s,t) S{x - Q{s,t))ds . 

This construction was extensively discussed in [18], where proofs were given in various cases. 
A key result is that the momentum map constructed this way is equivariant, which means it is 
also a Poisson map. This explains why the coordinates (Q,P) undergo Hamiltonian dynamics. 
There is also a right action (Qg = Qo g), whose momentum map corresponds to the canonical 
one-form P • dQ on T*Emb. 

The EPDiff equation has been applied in several contexts, ranging from turbulence model- 
ing [8] to the design of imaging techniques [211 I22j . Its CH form in one dimension has been 
widely studied because of the integrability of its singular solutions (peakons). The properties 
of geodesic flows on the diffeomorphisms regarded as a Lie group have played a central role 
in determining the behavior of the singular solutions of EPDiff. This suggests that we pursue 
further investigations of the emergence of singularities in the solutions of geodesic equations on 
Lie groups. 



1.3 The two-component Camassa-Holm system 

In recent years, the Camassa-Holm equation has been extended [6l[7l[23] in order to combine the 
integrabihty property with compressibihty, which introduces a pressure term in the equation for 
the fluid momentum. The resulting system (CH2) is a geodesic motion equation on Diff(DJ^, 
given as an Euler-Poincare equation on the semidirect product Lie algebra X (S) ^. In the general 
case, the Euler-Poincare equations are written on the dual of a semidirect product Lie algebra 
0(Dy as dH] 

whose components are 

d SL ,, SL 5L d SL ^ 5L 

-r- -rr = — adf -r— + -^ o a , -— -^ = — t, -^ 

dt oE, oE, da dt da da 

where the notation E,SL/6a stands for the (left) Lie algebra action of g on 1/*. 

The integrable CII2 equations are derived from the following variational principle on X(DJF 



/ L{u,p)dt = with L{u) = - iu{l-dl)udx + - I p^ dx. 



6 ' 
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Explicitly, the CII2 equations are 

Pt = - {pu)x , 

lit Uxxt -~ oUUx ~r ^"i^xUxx ~r UUxxx PPx • 

These equations describe geodesic motion with respect to the H^ metric in u and the L^ metric 
in p. Legendre transforming yields the metric Hamiltonian 



H{m, p) = — I m{l — d ) rridxH — I p dx . 



Extending to more general metrics (Green's functions) yields the Hamiltonian 

H{m, p) = - j j m{x) Gi {x — x) m{x') dx dx' + 7: p{x) G2 {x — x') p{x) dx dx' . 

This Hamiltonian yields the Lie-Poisson equations 

Pt = - {pu)x , 
mt = - unix - 2mux - pXx , 

in which u and A are defined by the variational derivatives, 

u = - — = Gi * m and A = -^ = G2 * p , 
dm dp 

and the symbol * denotes convolution. For example, we may write the Euler-Poincare equations 
for the H^ norms associated to the two length-scales ai and a2 (i.e., Gj = (1 — a^5^)~^) as: 

Xt - a\\xxt = -{u\- aluXxx)^ , 

ut - a\uxxt = -^uux + 2a\uxU^x + ct^uuxxx - Xx{\- alXxx) ■ 



The second equation may also be expressed in hydrodynamic form as 

ut + uux = —Px , with pressure p := Gi * iu + —u^ + -A — — A^, 

The integrable case CH2 is recovered from these equations for 02 = 0, although this case does 
not allow for singular delta-like solutions in the density variable (or depth variable, in the shallow 
water wave interpretation), since \ = p when 02 = 0. 

One may also extend the problem to higher dimensions by taking the Hamiltonian 

H{m,p) = - /"/"m(x)Gi(x-x')m(x')d"xd"x' + - /"/"/)(x) G2(x - x') /9(x') d"xd"x' (1.8) 

which represents the two-component extension of EPDiff and is denoted by EP(Diff(s)^). In 
three dimensions, the corresponding Lie-Poisson equations assume the form 

pt= -V ■ (pu) , 
lut = u X curl m — V(u • m) — m(div u) — pVA , 

where u = Gi * m and X = G2 * p- These expressions can also be written in a covariant form 
by using the Lie derivative £u with respect to the velocity vector field u G jC(M") of a one-form 
density 

m = m-dx®d"xeX*(M"), 

whose CO- vector components are m. The diamond operator (o) corresponding to the Lie 
derivative £u is defined by 



P,-£uX) :=(poA,u) , (1.9) 

where in the present case A € 1^* = ^ is a scalar function while p € V = Den is a density 
variable and ( • , • ) is L^ pairing in the corresponding spaces. In this notation, one writes the 
general case for p & V and A E F* as 

pt + £uP = 0, 
mt + £u'm = poX. (1-10) 

Upon recalling that the metrics Gi are Green's functions associated to the diff'erential operators 
Qi (that is, Qi • Gi = 5{x — x')), one finds the following Lagrangian 



L(u, A) = ^ /u Qiu d"x +^ I'xQ^X d"x 



whose Euler-Poincare equations recovers the EP(Diff(s)J^) equations p.lU|) . upon substituting 
the Legendre transforms m = Qi u and p = Q2 X. 

Specializing the operator Q2 to Q2 = A (that is, taking a2 — > 00 above) allows this system 
to be interpreted physically as an incompressible charged fluid. In this case, p is interpreted 
as charge density, rather than mass density, and A is the electrostatic interaction potential, 
satisfying A A = p for the Coulomb potential. The pure CII2 case {Q2 = 1, or a2 — > 0) 
corresponds to a delta-like interaction potential, which also appears in the integrable Benney 
system [3l[10]. The application of the full semidirect-product framework to the imaging method 
called metamorphosis was studied in [22]. The relation of the latter work with the current 
investigation will be traced further in the sections below. 



1.4 Plan of the paper 

This paper follows [18] in studying the singular solutions of geodesic flows on semidirect-product 
Lie groups by identifying them with momentum maps. The new features here are: (i) we explain 
their geometric nature in the context of the Kaluza-Klein formulation of a particle in a Yang- 
Mills field and (ii) we extend the applications of this approach to systems of equations governing 
geodesic flows on semidirect-product Lie groups with a non-Abelian gauge group. 

The next section first shows how this works when considering the Abelian gauge group of 
functions T . Momentum maps for left and right actions of the diffeomorphisms and the Lie 
group of gauge symmetry are derived that recover the corresponding conservation laws. The 
key observation is that the collective Hamiltonian for the singular solution momentum map is a 
Kaluza-Klein Hamiltonian, thereby recovering the conservation of the gauge charge. In section 
3 we show how this observation extends to the consideration of a non-Abelian gauge group G 
for particles carrying a spin-like variable. The introduction of a principal G-bundle becomes 
important in this context, since it best clarifies the reduction processes that are involved in 
the collective Hamiltonian system. Again, the collective Hamiltonian arises from a Kaluza- 
Klein formulation, after reduction by the gauge symmetry. We show that the Kaluza-Klein 
construction allows one to extend all the geometric features found in [18] for the EPDiff equation 
to the semidirect-product geodesic flows under consideration. 

The last section is devoted to illustrating how the continuum geodesic equations on semidirect 
products also arise naturally from a Kaluza-Klein formulation. In this formulation, the advection 
equation of the gauge charge density is recovered as the conservation of the momentum conjugate 
to a cyclic variable. The Kaluza-Klein approach also leads to a Kelvin circulation theorem in 
the continuum description. 

2 The singular solution momentum map for EP(DifF(S).^) 

Given the semidirect-product Lie-Poisson equations (jl.lOp . direct substitution shows that they 
allow the singular solutions 

N 

{m,p)=Y,l{Fi{s,t),w,{s))5{^-Q,{s,t))d''s, (2.1) 

where s is a coordinate on a submanifold S of M", exactly as in the case of EPDiff. If dim 5=1, 
then this corresponds to fluid variables supported on a filament, while dim 5 = 2 yields sheets 
of fluid density and momentum. The dynamics of (Qj, Pi,Wi) is given by 



dQiis,t) 
dt 



Y, j^jis',t)GiiQiis,t) - Qjis',t)) d'^s' . 



^^^= -^y"Pi(s,t)-P,(s',t)VQ.Gi(Qi(5,t)-Q,(s',t))dV 
j 

-^ Li(s)t/;,(s')VQ,G2(Q,(s,i)-Q,(s',i)) dV, 
j •' 

with dtWi{s) = 0, \/i. 

Recalling the geometric nature of the pulson solution of EPDiff and following the reasoning 
in [18], one interprets Qj as a smooth embedding in Emb(S', M") and Pj = Pj • dQj (no sum) as 
the canonical one-form on T*Emb(S', M") for the i-th pulson. In the case of EP(Diff(D.F), the 



weights Wi for i = 1,. . . ,N are considered as maps Wi : S ^ M*. That is, the weights Wi are 
distributions on S, so that Wi € Den(S), where Den := J^*. In particular, we consider the triple 

{Qi,Pi,Wi) e r*Emb(5,M") x Den(5) . 

and prove the following. 

Theorem 2.1 (Singular solution momentum map) The singular solutions h2.1\) of the semidirect- 
product Lie-Poisson equations I11.10\) are given by 

N 

(m, 






These expressions for (m, p) S X*(]R") (S) Den(IR"') identify a momentum map 



N 



J : X ('r*Emb(5, M") x Den(5)) -^ 3e*(]R") (S) Den(M") . 

Proof. For convenience, we fix label i and suppress the summations in our singular solution 
ansatz. In order to define a momentum map we first need to establish a Poisson structure on 
T*Emb X Den. Since the weights have no temporal evolution, it is reasonable to propose the 
canonical Poisson bracket on the new phase space, so that 

r^^.^^^ ^ ^ f f dF dG OF dG\ ,. 

{f.G)(Q,p..) :=Y.J [aQiwr - 3QJW,) ""■ 

Now, if /? = (/3i, /3o) G X(R") © T{W), then the pairing (J, (3) is naturally written as 

(j(Q, P,w),(3) = J (P(s) . /3i(Q(s)) + w{s) /?o(Q(s))) dh . (2.2) 

Consequently, one calculates 

'd(3j ^ dl3o\ 5F 



{F,(P,f3,{Q))} + {F,{w,(3oiQ))} =/(^-/3i(Q(^)) 



P+ w^] ■--] d^s. 



dQ dQj <5P 

(2.3) 



This may be written equivalently as 

{F,{3,p)}=Xp[F], 
in which the vector field Xp has (Q,P,i(;) components 

This vector field is identified with a Hamiltonian vector field corresponding to the Hamiltonian 

H = I (wis) (3o{Q{s)) + P(s)-/3i(Q(s)))d'=s = (j,/3). 

This Hamiltonian corresponds to compositions of cotangent lifts T*Diff generated by f3i with 
fiber translations 

T-d{wi3o) ■ (q^p) '■= {q,p-wdPo) 

generated by —wj3q (notice that w is independent of q, so that d(w;/?o) = wd/3o)- Thus Xp is 
an infinitesimal generator. ■ 

Remark 2.2 Theorem \2. 1\ proves that the singular solution is a momentum map deriving from 
the infinitesimal action o/ Diff (s) .F. However, it does not identify the corresponding global 
group action. This is explained in the next section below. 



2.1 A left group action for the singular solution momentum map 

The proof of Theorem 12.11 shows that the momentum map J m equation (j2.2p is obtained by 
the following left action of the semidirect-product group Diff (s) J- 



QW, PW, «;Wj = (77t(Q(°)), P(°) • r7?,-i(Q(°)) -d(u;/?o(Q(°))), u;^ 
= (7?,(QW), PW • Vr?,-i(QW) - V(u;/3o(Q(°))), u;^")) 



where (t/j ,/3o) € Diff ©JF. It is worth noticing that the order of operations in the composition 
is not relevant, since 



C^or^,,„ (q(o), P(0), ^(0)) = (^,,;3„)-i (q(o), p(0)^ 



w 



(0) 



so that (r\ ^,/3o) = ii]i~Po) ^ ^ Diff(s)jF. This can be easily seen from the following calcula- 
tion, where we take z/; = 1 for simplicity, 

{v,Po)-' o {rj,(3o) . (Q, P) = {v~\-Po) ■ (r/t (Q), P • Tr?,"! (Q) + d/3o(Q)) 

Consequently, exchanging the order of operations in the composition simply yields another 
element of Diff (s)^ and the arguments above are still valid. Such group operations are useful in 
deriving fluid descriptions from kinetic equations, where the so called plasma-to-fluid momentum 
m,ap determines the Hamiltonian structure of the fluid system [30] . 



Remark 2.3 (Geodesic motion on fiber-preserving transformations) The arguments 
above show that £^P(Diff (S).F) is a geodesic m,otion on the Lie group Sympjr(T*Emb) of fiber- 
preserving symplectomorphisms on the cotangent bundle T*Emb. It is a general result that 
any transformation by Symp^(T*Emb) can be realized as the composition of a fiber translation 
and a cotangent lift f^. On the other hand, such a transformation is always a canonical trans- 
formation characterized by a generating function that is linear in the canonical momentum. 
Therefore, EP(DiS(§)T) ~ EP Symp^ . Following this line of reasoning, one may also consider 
geodesic flows on the entire symplectic group Symp (T*Emb) . These flows apply to the moment 
hierarchy in kinetic theory and have been treated in 113\, \14^ . 

Another important property of the momentum map J in ()2.2p is its equivariance [26] , which 
guarantees that J is also a Poisson map. Below we show the infinitesimal equivariance of J, 
which is defined by the following relation 

X;3KJ,7>] = (J,ad^7) V/3,7GX(S)J-. 

Theorem 2.4 (Equivariance) The singular solution momentum map J in [K 

J(Q, P, w) =j{P{s, t),w{s)) 5(x - Q(s, t)) d^s 
is infinitesimally equivariant. 



10 



Proof. 






^Q^d^wyM.)^^,.^-7.-^)d«. 



j (jv{s,t)5{^-Ci{s,t))dh 
f ( fw{s)6{^-Q{s,t))d''s 



+ 

(J, ad/3 7) • 



(/3i • V)7i - (7i • V)/3i 
(/3i • V)7o - (7i • V)/3o 



d"x 
d"x 



2.2 Collective Hamiltonian 

From the expression of the vector field (12. 4p , we can immediately write the equations of motion 
for Q and P. Moreover, if we insert the singular solution momentum map into the Hamiltonian 
p^ we recover the following collective Hamiltonian Hn : x^^^ (T*Emb(S',M") x Den(5)) -^ M 

N 



Hn = ^^ /"/"Pi(s,i)-P,(s',t) Gi{Q,{s,t)-Qj{s',t)) d'^sd^s' 
hi '' '' 

N 
+ \iZ [[Ms)wJis')G2{Q^is,t)-Q,is',t)) d^'sdV. 



The corresponding canonical Hamiltonian equations of motion are 

dQi 
dt 

dPi 



kJ 



dt 



Y, / P.(^',i) Gi(Q,(s,t)-Q,(s',t)) dV, 

-Y^ J w,{s,t)wjis',t) ^G2{Qiis,t) - Qj{s',t)) dV. 



Remark 2.5 (Singular potential terms and integrable cases) In the limiting case G2 = 
6, we recover the integrable CH2 system corresponding to a geodesic motion with a delta-like 
potential. In terms of particle dynamics, this means that for a positive potential, two particles 
will bounce off immediately before colliding; while for a negative potential, they will proceed 
together, attached one to the other and they will never split apart. This kind of singular delta- 
like potential is also present in another integrable system, the Benney equation |3, IJOJ /. 

2.3 Kaluza-Klein formulation 

It is interesting that the collective Hamiltonian allows for a Kaluza-Klein formulation similar to 
the usual treatment of a particle in a magnetic field [26]. In order to see this, one first observes 



11 



that the equations of motion for {Qi,Pi,Wi) may be recovered on the Lagrangian side via the 
Legendre transform, 

This yields the fohowing Lagrangian L : x^^^ (rEmb(5,M") x T{S)) -^ R 

N 

2 



TV 

+ \t^ll 0'{s)e^{s')G''{Qiis,t)-Q,{s',t)) d'sdh' 

where G* with raised index is the inverse metric associated to Gi . (If Gi is given by a convolution 
kernel, then G* becomes a differential operator.) As our notation may suggest, we now enlarge 
our configuration space so that the Lagrangian Ljy becomes defined on 

TQkk := r(Emb(5,M") x J^{S)) = rEmb(5,M") x TT{S) 

and Qkk is called Kaluza-Klein configuration space. Now, since the coordinates 9i in the 
Lagrangian Lj^j = L(Q*, V*, ^*, 0*) are ignorable, its conjugate momenta wi will be constants 
of motion. This allows the collective Hamiltonian to be naturally written in the Kaluza-Klein 
formulation on T*Qkk- In this framework, it is well known [31] that the weights Wi are another 
type of conserved momentum map. 

Remark 2.6 (Physical interpretation) The physical system described by the Kaluza-Klein 
Hamiltonian H^ turns out to be related to the motion of electrical charges whose mutual interac- 
tion is given by the potential term in G2 ■ This relation is evident by noticing that for the case of a 
single particle, the Hamiltonian Hjy reduces to the Kaluza-Klein Hamiltonian Hi = P^/2-|-u;^/2 
of a free charge \26^ . In the multi-particle case, the momentum of a single charge is affected not 
only by the momenta of the remaining particles in the system, but also by their charges involved 
in the potential term. 

The absence of an external magnetic field in the system is reflected the absence of its vector 
potential A in the Hamiltonian H^. One may imagine inserting a non-zero magnetic field into 
the system by substituting Pi(s) — > Pi(s) — Wj(s)A(Qj(s)) into the Hamiltonian Hj^f. This is 
the usual process in electromagnetic (or Yang-Mills) theory. 

2.4 A right action momentum map and the Kelvin-Noether theorem 

As shown in the previous section, the singular solutions identify a momentum map which is 
determined by a left action of the group Diff(R") on Q G Emb(S', M"), i.e. r/Q = 77 o Q. 
However, as in the case of EPDiff, one may also construct a right action by Qry = Q o r/, which 
is defined through the group 015(5), rather than Diff(M"). 

In order to perform such a construction, we consider the Kaluza-Klein formulation from the 
previous section, so that the configuration space is now Qkk = Emb(5, M") x Den(5'). The 
group Diff (5) acts on Qkk from the right according to 

The cotangent lift of this action to T*Qkk yields the following. 
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Theorem 2.7 The map 

Js(Q, P, d, w) = P(s) • dQ(s) + w{s)de{s) (2.5) 

is a momentum m,ap 

Js : r*(Emb(5,M") x J^(5)) ^ X*{S) 

corresponding to the cotangent lift of the right action o/ DiflF(S') on Emb(S', M") x J-{S). This 
quantity is preserved by the flow generated by the Hamiltonian H^. 



Proof. Although we might proceed analogously to [18] by applying the formula for momentum 
maps arising from cotangent lifts, here we prefer to proceed by applying the general definition, 
since it exhibits the infinitesimal right action on each of the variables. Inserting the functional 

(js, /3i) = ^(P • V,Q + wV^e) ■ /3i(s)d^s 

into the Poisson bracket {F, (J,/3i)} yields 

{F,(Js,/3i)} = j[(3i{s) • VsQ • ^ + V, • (/3iP) • ^) dh 

+ I (/?i(s) . V,0 ^ + V, • {w /?i) ^) d'^s = ip, [F] 

where ^^3^ is the infinitesimal generator corresponding to cotangent lifts generated by /3i G X{S). 
The last part of the statement follows from the fact that the Hamiltonian H^ is invariant 
under the cotangent lift of the right action of Diff(S), which amounts to invariance of the integral 
over S under reparametrization using the change of variables formula. ■ 

Corollary 2.8 (Kelvin-Noether theorem) 

Exactly as in the case of EPDiff, the conservation of the momentum map in Theorem 2.1 
leads naturally to the Kelvin-Noether theorem. Indeed, upon recalling the evolution of w, i.e. 
ti; + Vs • (it; /3i) = 0, one checks that 

and therefore 

I ^ {w-\s) Pis) ■ dQ{s) + de) = j^j w-\s) P(.) • dQ(s) = , 

where "ft is any Lagrangian circuit and the first step is justified by Stokes ' theorem. 

Remark 2.9 (Dual pair structures) Upon recalling from \18^ that the term Ji := P • dQ 
in 112. 5\) is a momentum map Ji : T*Emb — > X*{S) for the right action of Diff(S), one may 
construct the same dual pair structure as in the geometric description of the EPDiff equation. 
Indeed, we may introduce the map 



Jsing(Q, P) = j^{s, t) 5(X - Q(S, t)) d'^S 



that is, the X* -component of the singular solution momentum map J in Theorem \2.1\ This is 
well known U8f to be a momentum map Jsing '■ r*Emb -^ X*(M") that generates the left leg of 



13 



the following dual pair picture 

T*Emb(5,] 




which is the standard dual pair picture associated to the EPDiff equation 1 18^ . Dual pair struc- 
tures have been used in 129^ to explore the geometric nature of Clehsch variables in fluid systems. 
In the present context the variables (Q,P) form the Clebsch representation associated to the dif- 
feom,orphism group DiflF(5). We refer to the original works flPl [^ for deeper discussions on 
the geometric nature of dual pairs. 

In the present case the right-action momentum map Js also takes into account the extra 
term J2 := wdO, which is associated with the space of scalar functions J'{S). This term is again 
a momentum map J2 : T*T -^ X*{S), which will be used later in this paper for the construction 
of another dual pair, associated to the J^* -component of the singular solution momentum map 
in Theorem\2.1[ 



3 Extension to anisotropic interactions 

This section extends the previous results on singular solutions to the case when the fluid motion 
depends on an extra degree of freedom, such as the fluid particle orientation. This occurs for 
example in the theory of liquid crystals [15] . A geometric fluid theory for such systems is already 
present in the literature regarding Yang-Mills charged fluids and quark-gluon plasmas |lll [T2] . 
This work formulates the equations for the fluid momentum m(x), mass density p(x) and charge 
density C(x), where the charge is considered as an extra degree of freedom of each fluid particle. 
(This is the colour charge in the case of chromo- hydrodynamics for quark-gluon plasmas.) These 
equations are written as 

dp , / 5H\ 

-- + div C ^ = ad^H C , 3.1 

ot \ omj jc 

dm ^ (5H \ / 5H\ ^6H /^^6H' 

IJF + ^ ■ ( to « ■") + r « feij ■" = - "^^ - ('^' ^«?, ,.,, 

where C takes values on the dual Lie algebra g*, whose corresponding coadjoint operation is 
denoted by ad*. Thus the charge variable C belongs to the space of g*-valued densities, which 
we denote by 

0*(M"):=Den(M")®0*, 

so that C £ fl*(M"). In what follows we shall use the elementary fact that the space fl*(M'^) is 
dual to fl(M") = J^ (g) g, and we shall use the same notation for q and g(]R"). The distinction 
should be clear from the different contexts. 

Now, the equations above are known to possess a Lie-Poisson Hamiltonian structure dual 
to the Lie algebra of the semidirect-product Lie group Diff(D(.F©g) [HI [l2]. One may also 
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consider geodesic Euler-Poincare equations on this semidirect-product Lie group. This problem 
has already been considered in [H] in terms of its singular solutions, although not in relation with 
momentum maps. For the sake of simplicity, we consider the semidirect product Diff ©g and 
denote the corresponding geodesic equations by EP(Diff ©g). (The commutative case g = .T^xM 
reduces to the case studied in Section O) 

In order to construct the EP(Diff (Dg) equations, one writes a purely quadratic Hamiltonian 



H{m,C) 



1 



m(x) • Gi(x - x') m(x') d"xd"x' 



+ 



C(x),G2(x-x')C(x'; 



d"xd"x' 



(3.2) 



which yields the geodesic equations in covariant form 

mt + £Gi*in in = Co{G2*C). 



(3.3) 



where the Lie derivative and diamond operations were introduced in Section [TJ In order to 
simplify the discussion, one can specialize to the case when particles have an orientation (or 
spin) in space and think of the charge density as the distribution of the local particle orientation 
in space, so that C € Den(]R"') 050 (3) ~ Den(R"') (8)M^ and the Lie bracket is given by the usual 
cross product. However, the following result applies in general. 



Theorem 3.1 The EP(Diff(Dg) equations admit singular solutions of the form 



N 

(m,C) = J] / (P,(s,t),/i,(s,t))<5(x-Q,(s,t)) d's; 



(3.4) 



associated with the momentum map 



N 



J : X (T*Emb(5,M") x g*(S)) ^ X*(]R") (S)0*(]R"). 

Proof. Again, we fix i for convenience and suppress it in the notation. Substitution of the 
solution ansatz (j3.4p into the EP(Diff(Dg) equations yields 



dt 



/3l(Q): 



dP 
'dt 



■P •Vq/3i-(m,Vq/3o>, 



'dt 



ad^o^ 



(3.5) 



In order to define a momentum map, we first need to establish a Poisson structure on r*Emb ^ g*. 
To this purpose we use the following Poisson bracket |lllll2j : 



{F,G}(Q,P,;.):=/(-.---.-Jd\. 



A*, 



6F_ 5G_ 
6n ' 6fi 



where (•, •) denotes the pairing on g* x g. Now, if /3 = {P^,(3q) € 
functional (J,/3) may be defined as 



d'^s, (3.6) 

g(]R"), then the 



(J(Q, P, /i), /?) = y ((m, /3o(Q)) + P • PAQ)) d^ 
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where (•, •) may now denote the pairing on either 0*x g or j£*(S)g*x X(S)0. (No confusion should 
arise from this notation.) At this point one calculates the Poisson bracket using equation ()3.6|) 
as, cf. equation ()2.3p . 



(3.7) 
Thus, one finds the Hamiltonian vector field {F, (J,/3)} = Xp[F], with (Q,P,/z) components 

The first two components of this vector field identify a Hamiltonian vector field on T*Emb 
corresponding to the Hamiltonian 



H = j ((//, /3o(Q)) + P-/3i(Q)jd'=s=(J,/3), 
which generates compositions of cotangent lifts T*Diff generated by /3^ with fiber translations 

'r-d(M,/3o> ■ (9'P) = {q^P- (M,d/3o)) 

generated by —{fi,f3o). The third component generates pure coadjoint motion of the charge 
variable // on g* according to 

,.(i) - \A* ^ ^ ,.(0) 

under the action of the Lie group G whose underlying Lie algebra is g := T^G. Thus, the 
three-component vector field X^ is an infinitesimal generator. ■ 

Remark 3.2 (Left action) Just as in the case o/ EP(Diff ©JT), we can write a similar left 
group action o/Diff ©g on T*Emb0g*. Indeed, by the arguments in the proof above, one sees 
that the momentum map J derives from the following left action 

(qW, PW, /.W) = (^, oQ(o), PW . r(r?,-i oQW) - d((^,/3o> o Q W), Ad;,p(_,,„) /.W) , 

where T]t = exp {t (3i) G Diff . The equivariance of J is proved via the same steps as we followed 
in the isotropic case. 

Remark 3.3 (Principal bundle structure of the configuration space) The proof above 
uses the fact that the space T*Emb x g* may be endowed with a Poisson structure. The ques- 
tion remains on how to identify a suitable configuration manifold corresponding to such a phase 
space. The answer can be found in the theory of gauged Lie-Poisson structures \24\ \25^ . such as 
the one we used in our proof. Indeed, it is well known how this structure arises from the reduced 
phase space of the configuration principal bundle B 

S = Emb(5,M") xG(5), 

where G{S) is the gauge group underlying g(5), consisting of G— valued maps g : S ^> G, so 
that TeG{S) = q{S). In what follows we shall use the same notation for G and G{S) and no 
confusion should arise from this choice. 
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It is also known that the phase space corresponding to such a configuration manifold is 
obtained by reduction of the cotangent bundle T*B. Thus our phase space can be obtained as 

r*Emb X 5* ~ T*B/G ~ (r*Emb x T*G)/G ~ r*Emb x T*G/G , 

where the last step is justified by the fact that G does not act on Emb. Therefore our geometric 
treatment follows from the principal bundle structure of the configuration space B. 

In the Abelian case studied in Section\^ one has G = T so that one may identify the gauge 
group with its Lie algebra. This identification is peculiar of Abelian gauge groups and cannot be 
performed in general. The next section shows how to treat the Kaluza-Klein formulation of the 
non- Abelian case. 

3.1 Kaluza-Klein collective Hamiltonian 

As in the isotropic case, we see again that the cohective Hamiltonian 

N 

2 



N » „ 

2 E/y (/^^(^)' G2iCl^is) - Qj{s')) f,,{s'))d''sd's 



^ , , , 

k I 



obtained by direct substitution of J(Q, P) in the EP(Diff ©g) Hamiltonian, allows for a Kaluza- 
Klein formulation. However, in this case the gauge group is not Abelian and we need to proceed 
more carefully. In the Kaluza-Klein picture of the motion of a colored particle in a Yang- Mills 
field, the particle motion is a geodesic on a principal G-bundle B. The metric on B is G-invariant 
and its geodesies are determined by the G-invariant quadratic Hamiltonian on T*B where the 
Poisson bracket is canonical. Specializing to our case yields 

TV 

Q^^ :=S= X (Emb x G) 
%=\ 

and since the second term in the Hamiltonian H^ is G-invariant by hypothesis, it may be lifted 
to T*Qkk as follows 



N 

2 



Hr,{Q\P^,g\pi) = IY1 [[Piis)-Pjis') Gi(QHs) -Q^(s'))d^'5d'=s' 



«j 



1 ^ 



«j 



where pi is the conjugate momentum of the group coordinate g^ € G, so that {g^,Pi) E T*G, 
and (•, •) is now the pairing between tangent and cotangent vectors on G. The momentum pi is 
conserved, since it is conjugate to the cyclic variable g*. Thus the Hamiltonian Hi\f is Klauza- 
Klein and thereby recovers the conservation of pi. Such a conservation law becomes coadjoint 
motion on the dual Lie algebra g*, such that 

Ai = ad|^^/5^^ f^i (no sum) , 

where /ij = g' ~^ pi Vi = l...A^ (no sum over i), exactly as happens for the motion of a rigid 
body, when G = SO{'i). As a consequence of the above arguments, it is clear how the dynamics 
(j3.5p of the singular solutions is Hamiltonian with respect to the Poisson bracket in (|3.6p , which 
is the sum of a canonical term and a Lie-Poisson term. 
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3.2 The right action momentum map and its imphcations 

One may also consider the right action through the group Diff(S'). Upon following the same 
procedure as for the abelian case, one finds the momentum map corresponding to the right 
action of Diff (S) 



Js(Q, P, g,p) = F{s) ■ dQ(s) + {^p{s),dg{s)) e r (5) (3.8) 

where the Kaluza-Klein phase space is now 

T*Qkk = T*Emh X T*G . 

This momentum map is again conserved because of the evident symmetry of the Hamiltonian 
Hjsf under relabelling s by a change of variables. 

Remark 3.4 (Kelvin-Noether theorem) 

Upon seeking a Kelvin-Noether theorem for the non- Abelian system, one recognises that this 
system does not provide any conserved density variable that could be used to construct the loop 
integral of a differential one form, as done for EP(DiS (^ T) in Corollary \ 2. 8[ For this purpose, 
it suffices to fix a weight w G Den(S') preserved by the flow to obtain the following circulation 
theorem 



d 
It 



i w-\s) (P(s) • dQ(s) + (p(s), dg{s))) = . 

As in Section 2, it is interesting to notice that the momentum map Js for relabelling symmetry 
by right action is determined by the sum Ji + J2 of two distinct momentum maps, one for Diff 
and the other for the gauge symmetry, 

Ji(Q,P) = P(s)-dQ(s) Ji :r*Emb^ je*(5), 

J2(5,P) = {p{s)A9{s)) J2 : T*G ^ T{S) . 

These momentum maps have the same target space, but different image spaces. Now, since 
the pairing {p{s), dg{s)) is invariant under the (left or right) G-action by cotangent lifts, it is 
possible to re-express it as 

{p{s), dg{s)) = {g~^p, g~^Vsg)ds=: {fi{s) , A{s)) ds = (^(s),^(s)) 

for a g-valued one form A{s) = A.{s)ds (i.e. a pure gauge connection). This result does not 
depend on the particular choice of left or right G-action, since the invariance property is not 
affected by this choice. Consequently, with the definitions fj, := g~^p and A := g~^dg one may 
rewrite J2 as 

J2(5,P) = {p{s),dg{s)) = {fi{s),A{s)) . 

Remark 3.5 The quantity A determines the following magnetic component of a Yang-Mills 
field 

B = d-^A = dA+ [A, A\ , 

which is localized on the embedded subspace S and is intrinsically generated by the moving charge 
fi. In fact, the connection A does not represent the Yang- Mills potential of an external force 
field (not present in this case), usually denoted by A = A(x) -dx, which is rather a g-valued one 
form over the physical space M". 

More formally, B is the curvature of the connection induced by A, which becomes important also 
in the geometric approach to the dynamics of complex fluids |3, \15f . In such an approach, the 
interpretation of A := g^^dg as a group one-cocycle plays a central role, as shown in |^. 
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The momentum map J2 can now be used to construct another dual pair, describing the 
geometry of the dynamics of the gauge charge /x. Indeed, it is well known [26] that the expression 
fi = g~^p =: 3R{g,p) is a momentum map Jr : T*G —^ g* associated to cotangent lifts of right 
translation. One may use this map to construct the following dual pair picture 



T*G{S) 




9*iS) X*{S) 

According to the general definition [32], a pair of momentum maps [)* < — P — > q* is called a 
dual pair if and only if Ker TJi and Ker TJ2 are symplectically orthogonal to one another. As 
explained in [12], a necessary condition for f)* < — P — > q* to be a dual pair is that each Lie 



group Gi associated to Jj acts transitively on the level sets of J^ with k ^ i. Now, Diff(5) acts 
transitively on the level sets of Jr = /i(s), because of the parameterization freedom. Moreover, 
the action of G{S) on level sets of J2 is transitive too, since it is given by cotangent lifts. Thus, 
similar arguments to those in [18] allow one to conclude that the above dual pair is properly 
defined. 

4 Kaluza-Klein equations for semi-direct products 

As we have seen from the previous sections, the Kaluza-Klein construction explains how collective 
motion on semidirect-product Lie groups arises under the singular solution momentum map. In 
this section, we extend the Kaluza-Klein formulation to continuum equations on a semidirect- 
product Lie group. The resulting equations apply to the method of metamorphosis in the 
problem of matching shapes using active templates in imaging science [22] . In this application, 
the zero level set of the momentum map for right action plays a crucial role. 

4.1 Case of Abelian gauge groups 

From the theory of semidirect-product reduction [27^ [28| I19j . one knows that the reduced La- 
grangian L : X(S).7^ ^ M is an invariant function on Diff ©.F. The reduction process may be 
performed in two ways. The first is to reduce according to the (right) action of the whole group 

T (Diff (s) J") / Diff (D J'^ ~ X (D j^ , 

while the second approach uses the fact that 

('t (Diff (s) ^) /.f) / Diff ~ X (D J^ , 

which is known as reduction by stages [5]. 

In order to construct the Kaluza-Klein Lagrangian for continuum semidirect-product motion, 
one must enlarge the semidirect-product structure to incorporate the cyclic variables. Indeed, 
one can consider the invariant reduced Lagrangian as a function 

Lkk -.XxTT^R, 

where the Lie algebra X acts on TT by the Lie derivative arising from the tangent lift of the 
(right) Diff-action on JF. In particular, we first consider the unreduced Lagrangian £(?/, ?), /, /) 
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on TDiff X TT. Since this is Diff- invariant, one may construct the Kaluza-Klein Lagrangian 
Lkk as 

^V,'n,fJ) = LKK{u,fr]~^,fr]~^) = Lkk{u,ip, X) , 

which then has exactly the same expression as the usual Lagrangian L on j£(D^, cf. |22j 

Lkk{u, v9, A) = - / u Qiu d"x + - / a Q2 a d"x = L(u, A) , 

although the function 99 is now considered as a cyclic variable, whose corresponding momentum 
(the density p) is preserved by the flow. 

From the arguments above we recognize that the Kaluza-Klein Lagrangian is related to the 
Lagrangian for EP(Diff (S)JF). Thus, it is natural to ask how this relation arises from a reduction 
process. One first observes that 

r(Diff (S) J')/Diff ~ (TDiff x rj^)/Diff ~ X x (TJ^/Diff) . 

The parenthesis in the last step indicates that functions in TJ^ are defined modulo the action of 
diffeomorphisms, so that any element (n, v) G TT/ Diff can be written as (/, /) 7]~^ \/i] € Diff. 
This application of the Kaluza-Klein reduction process to the unreduced Lagrangian C differs 
from ordinary semidirect-product reduction, which proceeds through the following steps 

(t (Diff (D T) /A I Diff ~ ('rDiff X TTIT\ / Diff ~ X X (.F/ Diff ) ~ X (D .^ . 

In conclusion, the two Lagrangians L and Lkk may be derived from the same configuration Lie 
group Diff (S)jr, by following different reduction processes. 

One may imagine that the same process may be followed upon replacing ^ by a generic 
vector space V, yielding the Lie algebra XCDTF. However, in such a case, the gauge Lie group 
needed to construct a Kaluza-Klein formulation is absent. 

4.2 Extension to non-abelian gauge groups 

All the considerations in the previous section may also be carried out for the case of anisotropic 
interactions, simply by replacing T by the non-Abelian gauge group G. The reduction process 
underlying the continuum dynamics on X(s)g proceeds as follows 

('t (Diff ©G) /J-)/ Diff ~ ('TDiff X TG/GVDiff ~X x g/Diff -XCDg. 

However, for G- invariant Lagrangians on X(S)0, we may interpret the dynamics as occurring on 
the space X x TG, so that the Kaluza-Klein Lagrangian 

Lkk : X X TG ^ M 

with (5,5) € TG and ?/ G Diff is written after right reduction by Diff as 

C{f]r]~^ , gr]~^ , gr]'^) = LKKiu,n,u) = 2 /" ^1" d"x + - {Q2i^, v) d"x, 

with definitions u := f]r]~^, n := (777"^ and v = grj^^- Now, since such a Lagrangian is also 
G-invariant, then one may write 

Lkk{u, n~^ n, n~^ u) = T(u, x) , 
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where n = gri~^ G G as before and L{\i, x) with x = t^~^ ^ is the Lagrangian on XCDg. Legendre 
transforming this Lagrangian produces the Hamiltonian H(m,C) in (j3.2p . This construction 
yields the conservation of the conjugate variable p = 5Lkk/5v along the flow of the group of 
diffeomorphisms, since n is an ignorable coordinate. The reduction process involved in such a 
system proceeds as follows 

T(Diff(DG)/Diff ~ (TDiff /Diff ) X (TG/Diff ) ~ X x (TG/Diff) , 

where TG is the group of tangent lifts of G, which is itself acted on by the diffeomorphisms. Thus, 
again the two Lagrangians L and Lkk niay be derived from the same unreduced Lagrangian C. 
Consequently, the geodesic motion on semidirect-product Lie groups of the kind Diff (S)G always 
possesses a Kaluza-Klein construction. 



4.3 Application to metamorphosis 

Lagrangian formulations on Diff(DG have been recently considered in [22], where the whole 
theory is extensively studied in the context of imaging science. The Euler-Poincare equations 
corresponding to a Lagrangian L(u, n, u) carrying the cyclic variable n = grj^^ € G are found to 
be 

5L _ 



d 



+ A 



dt 



5L 



6v 



du 



6u 



n 



(4.1) 



in which the last equation arises from the partial time derivative of the definition n = gri~^ G G. 
These equations imply that the Legendre-transformed variable p = 5L/5v is preserved by the 
flow, which does not occur in the general case, when L depends also on n. 

In the general case, one may obtain the dynamics directly from a constrained variational 
principle 5S = Q with 



S 



L{u,n,i^) + {p, — + £un 



dt, 



where the angle bracket denote L^ pairing. Stationary variations produce 
= 6S = 



pon,d\i) + {dp, — + £un 



+ 



SL 



+ 






dt. 



(4.2) 



in which the diamond operator (o) is defined via the natural generalization of (jl.9p and £u 
with superscript dagger denotes the L^ adjoint of the Lie derivative so that, in particular, 
{£uP, Sn) = {p, £u6n). Standard Euler-Poincare theory for the case that n is a scalar function 
and its dual p is a density then implies the following system after a brief calculation. 



li-'- 



(5L 



d_ 

dt 
d_ 
dt 



+ £u 



5L 

n = 



5L 
6n 



= 0, 
6L 



P, 



(4.3) 
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System (|4.3|) possesses an exchange symmetry between the variables (n, z^) € TG and their dual 
variables {6L/6n, 5L/5v) G TG* , and it satisfies the following proposition. 

Proposition 4.1 System ( [^.3| j is equivalent to system {■i.l\ j when bhjbn = 0. 

Proof. This proposition follows from a direction calculation using the chain rule for the 
diamond operation and substituting the last three equations in system ()4.3p . Namely, 

d \6L 6L 

When n is a scalar and p is a density, then — pou = {p, dv) and setting 5L/5n = recovers the 
first equation in the system (|4.ip . ■ 



Remark 4.2 The Legendre transformation of the constrained Lagrangian defines the Hamilto- 
nian 

H{in,p, n) = (m, u) + {p, v) - L(u, n, v) , (4.4) 



in terms of the fiber derivatives 



m := -—- ana p := -^ . (4.5) 

OVl ov 



The variational derivatives of the Hamiltonian are found by substituting into the Legendre trans- 
formation ^.4\ ) o,s 

/^H , \ /5H ^ \ /5H ^ 
6H = /^,6in) + {^,6p) + {^,6n 
\dm I \ dp I \ on 

u, (5m j + (v, 6p) - (—, 6n 

Hence, 6H/6ni = u, 6H/6p = v and the semidirect-product Lie-Poisson Hamiltonian equations 
corresponding to the system CTTl) are 



d ^ \ 6H 6H 

Ol + £sH/s.n)m = - — on + po — , 

§-^ + £sH/Sm)p=-^-^, (4.7) 

g-^ + £5H/Sm)n-—. 

Perhaps not unexpectedly, the first equation may also be written to agree with the motion equation 
in (4^ as 



\di^ -^^H/Sm j m = ( ^ + £5H/5m j (P <> ?^) ■ (4- 
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4.4 The Kelvin circulation theorem 

Proposition 14. II allows the Kelvin circulation theorem for these semidirect-product systems with 
5L/5n = to be expressed in two ways, upon introducing a conserved density variable p, 
satisfying 

On one hand, the first equation in the system ()4.ip and the general theory of dynamics on 
semidirect-product Lie groups [19] imply, with definitions (j4.6p . that 



-^/^ = -/i(p,d.>, (4.9) 

dt J^, p Li,p^ 

where 74 is a closed loop moving with the flow of the velocity vector field. On the other hand, 
the first equation in the system (j4.3p implies that the Kelvin circulation theorem may also be 
expressed as 

|£i(m + <p,dn))=0, (4.10) 

where the sum vntot := m + {p, dn) is the total momentum. This circulation theorem for total 
momentum is the natural extension to the continuum description of formula (j3.8p for preservation 
of the right-invariant momentum map. The two circulation laws (|4.9p and (j4.10p are shown to 
be equivalent in the following. 

Proposition 4.3 The two forms of the Kelvin circulation theorem in equations \4-^ o,nd \4.10() 
are equivalent. 



Proof. This statement will hold, provided 

— * - (p, dn) = d) - (p, di/) , 
dt L, p' /y, p ' 



which may be verified directly, by using the second two equations in the system (|4.ip as 



1/J^^''"> = 


£^ si"' <"•*•>) = £ ^"' (!<"•*'> ^^'■<^' "" 


= 


£p,"'((«^^"'''*')"'(^- '«"'''""")) 


= 


£p,"'((^''*^"""» 


= 


£^"'<''''"'> = £;^ <"•"">■ 



where one recalls that the exterior differential commutes with the Lie derivative and the last 
step follows from the equation for dn/dt. ■ 

This proposition extends the arguments in Section [3.21 to the non-Abelian case in the continuum 
fluid description. 



Remark 4.4 The zero level set of the total momentum, cf. equation ( tg.^j) . 

m + (p, d?i> = , (4.11) 
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is preserved by the first equation in the system j^. j[ j. The preservation of zero total momentum is 
a key step in the metamorphosis approach using active templates in imaging science, because the 
zero value is imposed by the requirement that an initial image would evolve to match a prescribed 
final image at a certain end point in time \22^ . 

The zero level set condition Ii4.11\ ) for total momentum imposes the relation 

m = - (p, dn) . (4.12) 

This is the equivariant momentum map obtained from the cotangent-lift of the right action of 
Diff on the gauge group G, defined by 

("^'">g*xg=(^^^'">rx0 = -<^' ^un)y.^. (4.13) 

Thus, the zero level set condition |^. J j[ j for total momentum is itself a momentum map. This 
particular momentum map also appears in the application of the classical Clebsch method of 
introducing canonical variables for fluid dynamics. See, e.g., flTi/ . 

5 Conclusions and open questions 

We have shown how continuum equations on a certain class of semidirect-product Lie groups 
aUow for singular solution momentum maps arising from the left action of diffeomorphisms on 
the G-bundle Emb(5', M") x G{S). On the other hand, the right action on Emb x G has been 
shown to yield another momentum map that recovers the Kelvin-Noether theorem. 

These results arose from the observation that the collective dynamics on Emb x G was 
generated by a Kaluza-Klein Hamiltonian, thereby recovering the conservation of a gauge charge 
from a cyclic coordinate in the gauge group G. 

The Kaluza-Klein construction for the collective motion was implemented in the continuum 
description by considering the semidirect product Diff (s) G as the product of the diffeomorphisms 
with a gauge group G. The Kaluza-Klein construction implies the Kelvin-Noether theorem. 

An important open question is whether the singular solutions (I2.ip and (13. 4p emerge from 
smooth initial conditions. For example, the Camassa-Holm equation shows the spontaneous 
emergence of singular solutions from any smooth initially confined velocity configuration [4] and 
it is reasonable to ask whether this feature is shared by the class of two-component Abelian and 
non-Abelian Camassa-Holm systems considered here. This question is being pursued elsewhere 

m- 

Another open question concerns more general semidirect products. In fact, the present dis- 
cussion has considered only semidirect products of the Diff group with G- valued scalar functions. 
However, in physical applications one may also find semidirect products of the form Diff (s)T, 
where T denotes tensor fields in physical space. The most important example is probably ideal 
magnetohydrodynamics, where T is the space of exact two forms (cf. e.g. J19j). The dynam- 
ics on such products differs substantially from the cases considered here and deserves further 
investigation. 
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